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From the summer cruise of the Mid-Atlantic Bight 
Experiment, conducted jointly by Naval Postgraduate School, 
University of Rhode Island, and Woods Hole Oceanographic 
Institution, a study of acoustic variability in the region 
of Mid-Atlantic Bight shelf break was conducted. The period 
of the experiment was from 19 July to 09 August 1996. The 
experiment consisted of a suite of acoustic and 
oceanographic sensors including three 400 Hz (100 Hz 
bandwidth) transceivers, one 224 Hz (16 Hz bandwidth) 
transceiver and two vertical line arrays (VLAs). This study 
involves the signal processing of data collected by a 
telemetry buoy, an analysis of the spatial and temporal 
coherence of the phones and plane-wave beams of the vertical 
array, and the tidal and seasonal variabilities of plane- 
wave arrivals at the vertical array. Results of the 
changes in arrival time of the beams, the horizontal 
displacement of the front, the changes in the speed of 
propagation of the wave, and the change in the water 


temperature are discussed. 
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I. INTRODUCTION 


A. OBJECTIVES OF THE MID-ATLANTIC BIGHT EXPERIMENT 

The end of the Cold War and the operations in Desert 
Storm have brought about new challenges in undersea warfare 
(USW) . The mission of USW is “battlespace dominance in 
support of Forward From The Sea” (Snyder, 1997). Navy 
systems must be technologically advanced and highly 
efficient in diverse littoral regions. To operate 
effectively, the Navy must perform research in these 
littoral regions. To meet this abi ec ey the Mid-Atlantic 
IB edele (MAB) experiment owaS conducted. The overall 
scientific goal of the MAB experiment was to quantify 
shelfbreak frontal variability and its coupling to adjacent 
slope water circulation and to determine the impact of this 
variability on sound propagation from the continental slope 
onto the shelf (Smith et al., 1996). The MAB experiment 
also focuses on determining temperature structure of the 
ocean environment via acoustic tomography. It will use the 
two aspects of tomography: the forward problem and the 
inverse problem. The forward problem emphasizes the forward 
propagation modeling and the generation of replicas to 
determine the limitations on how replicas affect tomographic 
methods. The inverse problem utilizes replicas to localize 
the environmental parameters. Furthermore, the MAB 


experiment has cultivated the study in the area of geo- 


acoustic inversions - extraction of bottom sound speed and 
density from data - and horizontal refraction (3-D 
effects). A secondary goal of the Mid-Atlantic Bight 
experiment is to obtain sufficient oceanographic and 
acoustical data to be able to contrast the wintertime and 
summertime variability of the shelfbreak front (Smith et 
al., 1996). The gain that the Navy gets from this study is 
fundamental insights into acoustic propagation in the slope- 
shelf eae aon region and application of these insights to 


target detection and localization. 
B. THESIS OBJECTIVES 


The objectives of this study are as follows. 
i Adapt signal processing algorithms to the real 
acoustical environment in order to process three days Ot 
real continuous vertical line array (VLA) data. 
Ze Quantify the variabilities of plane wave arrivals and 
compare these variabilities with expected variations in the 
environment. 
cr Analyze the spatial and temporal coherence of the 
phones and beams of the vertical array. 
a Determine whether the variabilities are due to 
nonlinear solitons passing through the region or larger 


scale frontal variability. 


II. OVERVIEW OF MID-ATLANTIC BIGHT EXPERIMENT 


The Mid-Atlantic Bight experiment (a combined 
oceanographic and acoustic experiment) was conducted jointly 
by the Naval Postgraduate School, the University of Rhode 
Island, and the Woods Hole Oceanographic Institution. The 
location of the experiment is the region of the continental 
shelf break south of Cape Cod. 

The Mid-Atlantic Bight experiment consisted of two 
cruises - a summer cruise (19 July to 9 August 1996) and 
Wintee w@eurse (09 February Co 27 Bebruary 1997). Theseummer 
and winter cruises were proposed to compare summertime and 
wintertime variation of the shelf break fronts. This study 
will focus on the summer cruise. The summer ee was 
aboard the R/V Endeavor from the University of Rhode Island. 
Major operations involved on the cruise were equipment 
deployment, operations of the Seasor, and hydrography. 
During the time of the experiment, the nature of the region 
was complex. The front was distorted by an eddy filament 
which was produced by a warm core ring absorbed by the Gulf 
Stream and the shelf water was moving soutwest (Pickart et 
ae. , Eo IG) 

The acoustical elements of this cruise were three 400 
HzeeSsources (100 Hz bandwidth). omen 224 12 souvee (16 shiz 
bandwidth), two 16 element vertical line arrays (VLA) or 


“tomography” arrays and one physical oceanographic array. 


The location of each acoustic element is listed in Table 1. 
The arrangement of the acoustical elements is illustrated in 
Fig. 1. In Fig. 1, the two "nenthernemost  tutametleseare the 
VLA's, the three triangles to the south are the three 400 
HZ sources. Tt is seen that the 400 Hz sources are on the 
seaward side transmitting to two VLA’sS on the shoreward 
Side. The concentration of this paper will be on the data 
received at the northwest (NW) VLA from the southwest (SW) 


400 Hz source at 300 meter depth. 


Latitude Longitude 
NW VLA 


ma ca al 
a. eal 


Table 1: Acoustic mooring locations for Mid-Atlantic Bight 
experiment. 
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Element VLA 





The Seasor is a towed mechanical fish provided by the 
Woods Hole Oceanographic Institution. Towed in a “yo-yo” 
manner behind the R/V Endeavor, it provided information on 
the water column properties during the experiment, including 
temperature and salinity as a function of depth. This data 
was later converted into sound speed profiles as shown in 


Fig. 2. This figure displays one data set taken along the 
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Figure 1: Arrangement of acoustical elements of the Mid- 
Atlantic Bight experiment. 
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Figure 2: Sound speed profile generated from Seasoar data 
collected during the MAB experiment (units of sound speed in 
m/s). 





track of interest in this study. The low sound speeds on 
the shelf are indicative of the cold water typical on 
northern latitude shelf regions. The strong shelf break 
fete OW) by from cold water to warm water, is also a 
characteristic feature of some shelf break environments. 
This large-scale front is expected to have a considerable 
influence on the acoustic propagation in the region and is 


one of the key features to be studied. 





IIIT. NUMERICAL METHODS AND SIGNAL PROCESSING 


A. MATCHED~-FILTER PROCESSING 

To obtain adequate signal-to-noise ratio (SNR) without 
losing resolution, modulated signals were transmitted 
extending the signal duration without reducing bandwidth. 
To correlate the received signal, a process known as 
matched-filtering is performed. The modulation of the 
Signal must be such that the autocorrelation of the received 
Signal will be much shorter than the transmitted signal. 
The matched-filter process is based on the assumption that 
the ocean sound channel can be modeled as a linear, time- 
invariant system. Based on this assumption, the received 
Signal is the sum of delayed replicas of the transmitted 


sagnal (Mum et alan, 1995) 


r(t)= Xk, s(t ~T,) (3.1) 


Were sake) MemomtEnesseceivedsscignal, k Gepthe amplitude er i) 
ray, Sit) is che transmitted Signal) and) wis Pea Sera! 


fame we Satine llea) « eatecny. The output of the matched-filter is 


given by 


p(t) =Jr(t)s *(t—T)dt, Ew) 


where p(t) is the output of the received signal after 


applying a matched-filter. The autocorrelation of Eq. (3.2) 


yields (Clay and Medwin, 1977) 


p(t) =| P(f)e°*df (3.3) 


where P(£) is™ehe product ors (ffamemwrif), the transmitred 
and received signal respectively, in the frequency domain. 
The matched-filter process, in general, correlates the 
received signal with a replica of the transmitted signal. 
The matched-filtering technique used to process the data 
from the Mid-Atlantic Bight experiment is a technique 
introduced by Birdsall and Metzger (1986). Birdsall and 
Metzger pioneered the idea of phase encoding a signal using 
binary finite-field algebraic ‘codes (also known as m- 


sequences). Using m-sequencing, a phase encoded signal can 


be written ac 


r(t) = Zcos(@,t + M(t)9,) (3.4) 


where r(t) is the phase encoded signal, Z is the amplitude 


of the signal, ®, is the carrier frequency, M(t) is the m- 


sequence to be encoded and 9, is the phase angle. The 


complex envelope of Eq. (3.4) 1s 


10 


S(t) = Be™% (3.5) 


where B is the amplitude of the signal. The autocorrelation 


is written as 
F(t) = | Be 0% Be MO og 7. (3.6) 


The acoustic signals of the MAB experiment consist of 
10 millisecond pulse transmissions. These transmissions were 
matched-filtered and coherently averaged over five minutes. 
They were then stored as transmission records at an interval 
of fifteen minutes. The fifteen-minute interval includes 
five minutes of transmitting and ten minutes of silence. An 
example of a transmission record is displayed in Fig. 3 asa 
depth-time plot. The phones are numbered vertically from 
top to bottom. The phone locations are from 31.5 to 84 
meterssan depth every 3.5 meters for a total of 16 phomee. 
The red region indicates the highest energy of the arrivals. 
B. PLANE-WAVE BEAMFORMING 

After an acoustic signal has been received, there are a. 
variety of techniques available to distinguish a particular 
feature of the propagation as well as reduce the influence 
of noise on the signal. Noise is usually generated near 
the surface and tends to occur in high modes or high angles 
of propagation (Jensen et al., 1994). Plane-wave 


beamforming is a process that involves decomposing a signal 
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into plane waves. This application 1s especially useful in 
underwater acoustics because the propagation of acoustic 
waves and the source/receiver locations are contained in the 
vertical structure of the acoustic field. Although the 
arrivals at the VLA will not simply be plane waves due to 
the variation in sound speed across the array, these 
variations are less than 2% and will be ignored here. 
Summing the outputs of spatially distributed sensors 
defines the concept of mean erm le A beamformer’s output 
is far more useful than an individual sensor for it allows 
one to listen preferentially to wavefronts propagating from 
one eke va over another. Therefore, a beamformer 
implements a spatial filter. For example, Defatta et al. 
(1988) describe beamforming using a line array of equally 
spaced elements. Let the position of each element be along 
the z-axis in an infinite, homogeneous medium with an 
undetermined number of remote sources. Then the output of an 


element at thesorigin ofthe coerdinate syoteniacemethea =m 
source 1s SCE): From the far-field plane-wave 


approximation, a source from direction 9, produces outputs 


Onabie seca 


fie) = sq{ e+ Sees) ). Ba) 


Cc 


where / is the array element index, d is the spacing of the 


m 


elements, @ is the angle of arrival from the m” source, and 


i 


c is the speed of propagation. To accomplish beamforming, 
one must apply weights and time delays to the individual 
element signals and then sum them coherently. The output of 


the beamformer in the direction 9, is then 


L ld\sin\@,,) — sin\@ 
b= 2% w,S,, ___lalsin(@,) — sin(@,)) (3.8) 


Cc 


where w are the element signal weights. 


To derive the directional response characteristics of 
the beamformer, consider a continuous, complex, plane-wave 


Signal propagating across the array, defined by 


SiC. ae (3.9) 


Because of the angle of arrival, 9, the source produces 


™m 


outputs of the form 


e,(t) — piPnt pilkmd sin(6,, ) 


(3.10) 


2 ,, ; th 
where k= 7 =—— is the wavenumber of the m_ source. Pie 
Cc 


m 


beamformer output is then 


M ; : z 
b (1) = ei Onl x we imknalsin(G5)-sin(On) W(6,) (3.1 1) 


14 


where 
L E 
WAAL) = 5 op, eptknteinen)-sini5m) (3.12) 


is the spatial Fourier transform of the array weights 


(Defatta et al., 1988). 


For all of the beamforming shown in this thesis, a 
beamforming program written in matlab was used. The program 
was tested on model data of Known characteristics before 
being used on real data. Comparisons of real beamformed 


data and model predictions were also made. 


soe Predicted Arrival Structures 


To predict the arrival structure of the acoustical 
beams, the Monterey-Miami Parabolic Equation (MMPE) model 
was used. This model, formerly known as the University of 
Miami Parabolic Equation (UMPE) model (Smith and Tappert, 
1995), uses the parabolic approximation method with a split- 
step Fourier algorithm (Hardin and Tappert, 1973) to solve 
the acoustic wave equation in cylindrical coordinates. 
Selected outputs generated by the model are displayed in 
Figs. 4 through 6. Three different environments were chosen 
to contrast the arrival structure predicted by the model. 
The first environment assumes an isospeed water column with 
c=1500 m/s and a constant depth of 100 m. The source depth 
for this case is 50 m. The second environment also assumes 


an isospeed water column with c=1500 m/s but uses bathymetry 
1.5 


extracted from the MAB region between the SW 400 Hz source 
and NW VLA. The third environment uses both bathymetry and 
sound speed data extracted from the MAB region. For the 
second and third environments, the source depth is 300 m. 
In all cases, the source is modeled as a point source with a 
center frequency of 400 Hz and 100 Hz bandwidth. Also, all 
three environments assume a homogeneous fluid bottom type 
with sound speed, density, and attenuation of 1600 m/s, 1.2 
g/cm, and 0.1 dB/m/kHz*, respectively. Figures 4 through 6 
show the arrival structures of the beams as a function of 
arrival time and arrival angle. In these figures, it is 
seen that changes to the environment affect Ehe vertical 
arrival structure at the array. The difference in arrival 
angle of primary energy 1s seen in Figs. 4 and 5 due to the 
upslope bathymetry input into the model. Specifically, 


energy arriving at grazing is reduced and the largest 


arrivals are observed near 10°. Next, when the MAB sound 


speed profile is included, the simple symmetry of arrivals 


in the previous cases is removed, as seen in Fig. 6. Most 


of the energy arriving between -10° and 10° is observed to 
occur within about them fiisetOeden Se Between +/- 10° and 


15°, the arrivals rapidly stretch over 0.5 s. This is 


typical of ducted profiles where propagation angles away 
from grazing encounter faster sound speeds which counter the 
path lengths. Beyond some angle, however, the increased 


path length dominates the travel time and arrival times 
16 
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Figure 4: Arrival angle-arrival time display of a simple 
i1sospeed, isobath environment. 
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arrival time display of isospeed 


Arrival angle 
environment with the bathymetry of the MAB region. 
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of the energy arriving between -10° and 10° is observed to 
Seeut Weendrn =aoout the first 0.1 s. Between +/- 10° and 


15°, the arrivals rapidly stretch over 0.5 s. AP clit een es 


typical of ducted profiles where propagation angles away 
from grazing encounter faster sound speeds which counter the 
path lengths. Beyond some angle, however, the increased 
path length dominates the travel time and arrival times 
increase rapidly. 

2. Observed Arrival Structures 

During the MAB experiment, approximately 90 gigabytes 
of data was acquired and stored on a telemetry buoy. This 
data was transferred from the telemetry buoy to 8-millimeter 
tapes. There were nine disks (01-09) of data available for 
processing. Of the nine disks, however, there are three 
gaps in the data due to equipment failures (Heydt, 1996). 
The longest period without gaps in the data is on disks 07 
through 9. These disks contain 3 days of data. aie) 
beamform the data, the depth-time matrix was read into the 
beamformer algorithm. A Hamming filter was used to reduce 
sidelobes of One Signal. The depth-time matrix was then 
zero padded (increase the number of points in the matrix to 
a value of 92) which prepared the matrisew tor @an Fria A 


series of FFT’s were performed to transform from time to 


frequency domain, depth to vertical wavenumber. space, 
vertical wavenumber space to angle, and finally from 
frequency back to the time domain. The end result is an 


Zu 


1aemeLcaleatO Figs i. The duration of the arrival is 
approximately 0.75 seconds, Similar to the prediction. The 
signals peak with early arrivals and drop off sharply to 


negligible amounts. 
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Arrival angle-arrival time display of single SW 


Figure 7: 
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to NW transmission that has been beamformed. 





IV. DATA ANALYSIS 


A. SPATIAL COHERENCE 


To analyze the spatial coherence of the phones, a 
Single record of the VLA data was examined. The peak 
correlation of each phone’s data with the bottom phone was 
computed and normalized such that the autocorrelation of the 
bottom phone was identically one. This calculation was then 
repeated relative to the top phone. Figures 8 and 9 show 
the phone to phone peak correlation of phones from bottom to 


top and top to bottom, respectively. These results show 


that the signal correlation drops to ~0.5 over a few phones, 


at ~10 meters in depth. This suggests different multipath 
arrivals are interacting across the length of the array at 
slightly different times. 

To observe the correlation of the beams, a single 
transmission was beamformed. After beamforming, the center 
beam was extracted from the angle-time matrix and the peak 
correlation with higher angle beams was computed, normalized 
such that the autocorrelation of the center beam with itself 


1s identically one. Similarly, the 10° and 20° beams were 


extracted and cross-correlated with higher angle beams. 
Figures 10 through 12 show these beam to beam correlations. 


In each case, the correlation begin to rapidly decrease 
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Phone To Phone Peak Correlation 


Phones (Bottom To Top 





Figure 8: Peak cross-correlation of the phones from bottom 
to top phone. 
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Phone To Phone Peak Correlation 


10 


Phones (Top To Bottom) 





Figure 9: Peak cross-correlation of the phones from top to 
bottom phone. 


Zp 


above ~27° and drops to a minimum for the beams above ~32”. 
This is due to a lack of arrival energy beyond this angle 


and is consistent with a local critical angle between 25° 


and 30°. With a water column sound speed at the bottom of 


1480 m/s, this suggests a local sediment sound speed between 
1630 m/s and 1700 m/s. This appears to be a bit large and 
is possibly due to some side-lobe leakage into angles which 
do not have any real incident energy. However, direct 


measurements of geophysical parameters at a nearby site do 


reveal a sound speed at the top of the sediment of ~1600 m/s 


with a significant gradient of ~20 m/s/m (Apel et al., 
997). 

The highest correlation levels are found in Fig. 10, 
the cross-correlation of the central beam with higher 
angles. Because of the “BetouE: BANE of the propagation, 
no significant arrivals are expected near grazing. This 
appears to be confirmed in the beamformed data displayed in 
Fig. 7 which indicates dominant arrivals at non-grazing 
angles. The high correlation values in Fig. 10 are then 
presumably due to sidelobe leakage from dominant arrivals at 
higher angles. Figure 11 shows a relatively rapid 


decorrelation between the 10° beam and higher angles 


consistent with a significant, unique arrival near 10°. 


Figure 12 shows a wider correlation peak. However, there 


does not appear to be a significant amount of energy 
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Center Beam To Beam Correlation 





Figure 10: Peak cross-correlation of center beam with beams 
ie to 40" . 
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Ten Degree Beam To Beam Correlation 





Figure 11: Peak cross-correlation of 10° beam with beams up 
Eo r0-. 
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Twenty Degree Beam To Beam Correlation 
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Figure 12: Peak cross-correlation of 20° beam with beams up 
Bowl) = 
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arriving at 20°, and so the high correlation over ~5° may 


also be an artifact of sidelobe leakage. 


B. TEMPORAL COHERENCE 

Temporal coherence of the bottom phones was completed 
by reading in 3 hours of data (12 transmission records and 
correlating the center/bottom phone of record 1 with itself 


and with the center/bottom phone of the other 11 records. 
The same was done for beams at angles 0°, 10°, and 20°. In 


all cases, the peak correlation was extracted by shifting 
the absolute time of each record. Figure 13 shows the 
temporal coherence of the center and bottom phones. The 


correlation 1S seen to drop rapidly from the first ping to 
the second and then maintain a relatively low value of ~0.5. 


This indicates that the correlation time of the acoustic 
Eransmission is less than our 15 minute sampling. Figure i 
shows the ping to ping correlation of the beams. Again, 
these plots show a rapid decorrelation of the signal from 
pang Tomping. In order to investigate the actual temporal 
coherence, it would be necessary to return to the original 


data before performing the 5-minute averaging. 
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Center Phone Ping To Ping Correlation 


Bottom Phone Ping To Ping Correlation 





Figure 13: Ping to ping peak correlation of the center phone 
and bottom phone over a 3-hour period. 
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Center Beam Ping To Ping Correlation 


Ten Degres Beam Ping To Ping Correlation 
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Figure 14: From top to bottom, “ping ste ping) peak cree 


correlation of the center, 10°, and 20° beams over a 3-hour 
period. 
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oe TIDAL AND FRONTAL TRAVEL TIME INFLUENCES 

In Figs. 15 and 16 the arrival structure of the center 
phone is plotted over a three-hour period (sampled every 15 
Minutes) and a three-day period (sampled every hour), 
respectively. Similar plots for other phones and beams have 
been examined and exhibit the same general characteristics 
of arrival time versus geotime, and so will not be 
explicitly presented here. Consistent with the rapid 
decorrelation computed in the previous section, Fig. 15 
shows the large ping-to-ping variation during a three-hour 
period. Despite these variations in the fine-scale 
structure over short times, a general trend in the early, 
peak arrivals can be observed over longer times, as 
exhibited in Fig. 16. Over the course of three days, the 


earliest arrival of each record changes from an arrival time 


Of ~2#775 euto -27.825 s, a difference of AT)=0.050 SmmErn 
addition, a periodic tidal response of ~12 hours is also 


observed with peak travel time fluctuations of At ~0.030 s. 


Let us consider the tidal response ye Physically, 
the primary influence of the tide on the shelf break front 
is to push the cold water “tongue” of the front water back 
and forth over the shelf break. As seen in Fig. 2, this 
protrusion of cold shelf water occupies roughly 50 m in 
depth or between % and ¥% of the water column beyond the 


ehnelt break. In Order to €CStimatevtne Meurzenta 
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display of arrival structure 


for center phone over a period of 3 hours. 
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Figure 16: Geotime-arrival time display of arrival structure 
for center phone over a period of 3 days. 
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displacement of the front, let us assume a Simplistic model 
Mim—vGitichina wsect tom Of wthe cold water front of horizontal 
extent R, is replaced by warmer water. Furthermore, we will 


assume the early arrivals are due to energy propagating at 
~10° from grazing, and we will neglect changes in the 


propagation path. Due to the near grazing propagation, the 


acoustic path length through this segment is of order % R 


£ 


(che % due to the limited vertical extent). int the “seund 


speed over this path changes from C,.,, = 1475 m/s to C,_ 


=1500 m/s, then the travel time will be reduced by an amount 


Naar a 
4 : Cspeif C 


warm 





). (4.1) 


Given At ~0.03 s for the tidal response, the horizontal 


displacement of the front is then R, ~10 km. This value is 


consistent with other oceanographic observations (Chiu, 


LC) The change in peak arrival times over the 3-day 
period examined, AT ~0.05 s, was due to a cold shelf water 


intrusion into the operational area during the experiment 


Ghai, ©1997). Again, let us assume a simplistic model in 
which all of the water on the shelf, covering the final ~10 


km of propagation, was affected by a reduction in sound 


speed f£rom’C..,, =1485 m/s tose 2, = 14es mee Dee Sornee 


this affects the entire path over the shelf, and since the 
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propagation is still nearly horizontal, we can estimate the 
change in travel time from 
l R hey 


ft 
DN ee —— )= 
fae C shelf C shelf Cor 








AGA: (4.2) 


imc AG. -lO0om/s whick® corresponds to a nearly 2°C change 


in the water temperature. This appears to be about twice 
what is expected, and so it may be speculated that the cold 
water intrusion significantly perturbed the propagation 
path, increasing both the propagation angles and the total 
path length. Note that this would be consistent with the 
observed decrease in arrival amplitudes since higher 
propagation angles would also generate higher bottom 


abienuation. 


or, 





V. CONCLUSION 


The Mid-Atlantic Bight experiment was conducted to give 
the U. S. Navy a better insight into shallow water acoustic 
propagation near shelf break regions. There was 90 
gigabytes of data collected during the experiment. Thass 
study involved the signal processing of a small portion of 
this data using the technique of plane-wave beamforming. 
The analysis of eae beams at the vertical array revealed the 
following results. 
de The phone to phone peak correlations show that the 


Signal correlation drops to ~0.5 over a few phones, at ~10 


meters in depth. This suggests different multipath arrivals 
are interacting across the length of the array at slightly 
different times. 

2s Beam to beam correlations show a rapid decrease for the 
beams above ~27°. thas 1s due to a lagheof arama aecnera, 
beyond this angle and is consistent with a local chai eo -Werul 


angle of between 25° and 30°. With a water column sound 


speed at the bottom of -~1480 m/s, this suggests a local . 


sediment sound speed of between 1630 m/s and 1700 m/s. The 
highest correlation levels are found in the  cross- 
correlation of the central beam with higher angles; however, 
because of the upslope nature of the propagation, no 


Significant arrivals are expected near grazing. 


Be, 


Sx The temporal coherence of the center and bottom phones 


is seen to drop rapidly from the first ping to the second 
and then maintain a relatively low value of ~0.5. Pais 


indicates that the correlation time of the acoustic 
transmission is less than our 15-minute sampling. Ping-to- 
ping correlations of the beams show a rapid decorrelation of 
the signal. In order to investigate the actual temporal 
coherence, it will be necessary to return to the original 
data before performing the 5-minute averaging. 

a Consistent with the rapid decorrelation of temporal 
coherence, there is a large ping-to-ping variation absorbed 
during a three-hour period. Despite these variations in the 
fine-scale structure over short times, a general trend in 
the early, peak arrivals can be observed over longer times. 


Over the course of three days, the earliest arrival of each 
record changes from an arrival time of ~ 27.775 s to = 1 oes 
s, a difference of AT ~0.050 s. In addition, a periodic 
tidal response of ~12 hours is also observed with peak 
Cravel time fluctuations of Aty_070J0n>. 

cog Given At ~0.03 s for the tidal response, the horizontal 


displacement of the front is then R, ~10 km, which is 
consistent with other observations. The change in the peak 
arrival times over the 3-day peroid is due to a cold shelf 
water intrusion into the operational area of the experiment. 


A simple analysis suggests AC,,.,, ~10 m/s which corresponds 
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to a nearly 2°C change in the water temperature. sUakks 


appears to be about twice what 1S expected, and so it may be 
speculated that the cold water intrusion significantly 
perturbed the propagation path, increasing both the 
propagation angles and the total path length. This would 
also be consistent with the observed decrease in arrival 
amplitudes since higher propagation angles would also 


generate higher bottom attenuation. 
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